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The origin of molecular hydrogen emission in cooling-flow filaments* 
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ABSTRACT 

The optical filaments found in many cooling flows in galaxy clusters consist of low den- 
sity (~ 10 3 cm~ 3 ) cool (~ 10 3 K) gas surrounded by significant amounts of cosmic -ray and 
magnetic-field energy. Their spectra show anomalously strong low-ionization and molecular 
emission lines when compared with galactic molecular clouds exposed to ionizing radiation 
such as the Orion complex. Previous studies have shown that the spectra cannot be produced 
by O-star photoionization. Here we calculate the physical conditions in dusty gas that is well 
shielded from external sources of ionizing photons and is energized either by cosmic rays 
or dissipative MHD waves. Strong molecular hydrogen lines, with relative intensities similar 
to those observed, are produced. Selection effects introduced by the microphysics produce 
a correlation between the H2 line upper level energy and the population temperature. These 
selection effects allow a purely collisional gas to produce H2 emission that masquerades as 
starlight-pumped H2 but with intensities that are far stronger. This physics may find applica- 
tion to any environment where a broad range of gas densities or heating rates occur. 

Key words: galaxies: clusters: general - galaxies: clusters: individual: NGC 1275 - galaxies: 
clusters: individual: NGC 4696 - intergalactic medium - infrared: galaxies 



1 INTRODUCTION 

The origin of the optical filaments often found near central re- 
gions of massive cl usters of galaxies remains mysterious (see 
. I Johnstone et allr2007l for a discussion of the extensive literature). 
The presence of dust and blue light, particularly in the centers of 
these objects, suggests that star formation has occurred, at least in 
some regions. Their spectra are quite different from H II regions, 
which we take as representative of emission from gas near newly- 
formed stars. Low-ionization optical forbidden lines and near-IR 
H2 lines are very strong relative to hydrogen recombi nation lines. 
No o ne model can account for these observations ( Johnsto ne et all 
l2007h . 

The combination of Spitzer and ground-based observations 
allows H2 lines from a wide r ange of excitation energies to 
be studied ( Johnston e et all 120071) . The resulting H2 population- 
excitation diagram shows that the level population temperatures 
correlate with upper-level energies. The high-excitation lines indi- 
cate r pop ~ 2000 K, substantially hotter than lines from low-lying 
levels, with T DOD ~ 300 K. A similar effect occurs in PDRs near 
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gal actic H II regions du e to starlight photoexcitation of higher lev- 
els ( Ivan Dishoeckl2004l) . The lines produced by this process are far 
weaker than is seen in these filaments. Generally, radiation from 
early-type stars cannot produce extremely strong H2 lines because 
the luminosity in the lines is, at most, a small fraction o f the stellar 
luminosity in the 1200 - 912 A range (Sternberg 2WSi). Other en- 
ergy so urces are needed. Shocks due to stellar outflows are one pos - 
sibility dFernandes et alll997l ; |jaffe et all200lMWilman et all2002h . 

At least two sources of energy exist that are unique to the 
cooling-flow environment. The filaments are surrounded by a mix 
of hot thermal particles and sy nchrotron-emitting relativistic par- 
ticles dSanders & FabiarJl2007l) . The latter are referred to as cos- 
mic rays here. Additionally, Faraday rotation measures and the 
long linear geometries suggest that m agnetic fields are important 
dGuidetti et alll2007l : lTavlor et alll2007l) . 

Here we investigate whether non-radiative heating, produced 
by cosmic rays or dissipative MHD waves, could produce the ob- 
served emission. Both would deposit energy into otherwise well- 
shielded molecular gas. We consider the cores of the filaments, 
ignoring complications such as a possible ionized sheath that 
would be produced by in situ star formation or mixing layers 
dCrawford & Fabian|[l992T) . We show that gas energized by non- 
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Table 1. H2 line intensities in a filament in the NGC1275 cooling-flow and 
the Orion Bar 



H2 Line 


I / I(Pa) flow 


I / I(Pa) Orion 


0-0 S(l) 17.03 \xm 


0.65 


lxlO -2 


0-0 S(2) 12.28 \xm 


0.30 


9xl0~ 3 


1-0S(1) 2.121 yun 


0.70 


6xl0~ 3 



radiative heating processes can produce an H2 spectrum that is very 
similar to those that are observed. 

2 CONDITIONS IN SHIELDED CORES 

Table Q] compares H2 line intensities relative to Pa in the Perseus 
cooling flow with those in the Ori on complex. Intensiti es for the fil- 
aments are taken from Table 8 of (Johnstone et al 2007) while those 
for the Orion Bar are from Sell gren etall ~Il990) and Allers et al 
(2005). The Orion Bar is an H + / H° / H2 interface viewed roughly 
edge on. Its successive layers of ionization are displaced on the sky, 
as shown in Section 8.5 of Osterbrock & Ferland (2006, hereafter 
AGN3). The separate regions are not spatially resolved in the extra- 
galactic case, so the intensities across the Bar were co-added to cre- 
ate the values listed in Table[T] This procedure is not more accurate 
than a factor of two. Even so, the differences in the line intensities 
are dramatic. The H2 lines are far stronger in the filaments than is 
found in this prototypical H II region. This suggests that additional 
energy sources are active in the cooling flow filaments. 

Could cosmic rays or extra heating, perhaps due to dissipative 
MHD waves, produce the strong molecular emission? The environ- 
ment inside a real filament is likely to have a mix of particle den- 
sities and be exposed to a range of cosmic rays, dissipative MHD 
waves, and possibly starlight. We consider only the cosmic ray and 
wave heating cases here. We do not address questions such as the 
transport of cosmic rays into possibly magnetically confined fila- 
ments, or the process by which kinetic energy in MHD waves might 
be converted into thermal energy. Rather we follow the effects of 
these energy sources on the microphysics of the gas and predict the 
resulting H2 line intensities. 

We use the develop ment version of the s pectral synthesis code 
Cloudy, last described by Ferland et al ( 1998). AGN3 discuss much 
of the physics of the environment with extensive reviews of plasma 
simulation codes. Our tr eatment of the H2 molecule, the focus of 
this paper, is described in lShaw"et a3( l2005l) . 

Our treatment of cos mic rays is des cribed in 

iFerland & Mushotzkvl h984l) , I Abel et all J2005I) , IShaw et all 
d2006h and I Shaw et all d2007h . In molecular gas a primary cosmic 
ray will produce a shower of secondary electrons which ionize 
and excite atoms and molecules. Relatively little of the energy 
goes into heating the gas (AGN3). Conversely, if the gas is highly 
ionized then nearly all the cosmic-rays energy will be thermalized 
and little ionization or excitation occurs. 

We also consider an "extra heating" case where thermal en- 
ergy, perhaps produced by dissipative MHD waves, is injected. This 
energy, ultimately due to the kinetic energy of a filament, is added 
to the thermalized kinetic energy of the particles in the material. 
We assume that the combination of plasma stream instabilities and 
particle-particle collisional coupling is effective enough in main- 
taining a Maxwellian distribution for the particle distribution func- 
tions that the detail of the energy injection process is unimportant. 
Ion-neutral drift, where charged particles are coupled to the field 



and m ove relative to neutral partic les, might be one source of such 
heat JLoewenstein & Fabian 1990). The effects are quite different 
from cosmic rays since the heating has no corresponding direct ion- 
ization or internal excitation. Ionization and excitation only occurs 
when the kinetic temperature becomes high enough to overcome 
energetic thresholds. Comparing the two cases, at a given kinetic 
temperature gas energized by cosmic rays will be more highly ion- 
ized and excited than the "extra-heating" case. 

We parameterize the cosmic ray case by the cosmic-ray den- 
sity relative to the galactic background value. Background cosmic 
rays produce an H° ion ization rate of 2.5 x 10~ 17 s -1 in gas with 
a low electron fraction IWilliams et all 19981) . The range of relative 
cosmi c ray rates includes the values measured by Sand ers & Fabianl 
d2007l) and is summarized in Section 3 below. In the "extra heating" 
case, energy is added to the thermal kinetic energy of the particles 
at a rate specified in erg cm -3 s -1 . It is not now possible to unam- 
biguously convert_ajnMSjTOdjriagnetic field into a wave dissipa- 
tion rate JHeiles & Crutchenl2005l) so we only show the range of 
rates with observable effects. 

H2 emission is likely to be produced in shielded cores where 
any ionizing radiation produced by either star formation or cooling- 
flow emission has been extinguished by surrounding gas. The 
z = metagalactic background is the radiation field incident upon 
the cloud. This includ es the CMB and the 2005 version of the 
Haard t & Madaul dl996h background with both starburst and quasar 
continua. This continuum is attenuated by a dusty gas with a total 
hydrogen column density of 10 21 cm -2 . The extinguished contin- 
uum has little effect on the simulations presented below but is in- 
cluded for completeness. For the "extra-heating" case we also in- 
clude cosmic rays with the galactic background value. 

There is no direct measurement of the gas metallicity or dust to 
gas ratio in the filaments. For simplicity we assu me galactic inter - 
stellar medium gas-phase abundances and dust ( Abel et "a1 l200l) . 
Many of the results presented below do not depend critically on 
these assumptions. Molecular hydrogen forms on grain surfaces in 
dusty environments. The rate of formation depends on the grain size 
distribution, temperature, composition and dust-to-gas ratio, all of 
which are unknow n. We adop t the galactic ISM catalysis rate of 
3 x 10 -17 cm 3 s -1 Jjurall 19751) . 

We are interested in molecular regions that must be well 
shielded from ionizing radiation. The absorbing column imposed 
on the metagalactic radiation field ensures that little light shortward 
of 912A penetrates. We further assume "Case B" (AGN3), that res- 
onance lines of H I and H2 have large optical depths because of this 
large column density. The result is the Lyman lines do not escape 
and continuum fluorescent pumping of H I and H2 resonance lines 
is not important. Absorption of UV radiation by electronic transi- 
tions is the main H2 destruction process in star-forming regions but 
will be unimportant in the calculations presented below because of 
the shielding provided by the large column density. 

With these assumptions the free parameters are the hydrogen 
density «h (cm -3 ) and the cosmic ray or extra heating rates. Pairs 
of figures, giving the cosmic ray case as the upper panel and the 
"extra-heating" case as the lower one, will be shown below. The 
hydrogen density is the x-axis in all cases. 

FigureQjshows the log of the gas kinetic temperature as a func- 
tion of these parameters. The temperature is the result of heating 
and cooling processes. In the cosmic ray case the heating per unit 
volume is proportional to rcR/m, where rcR is the primary cosmic 
ray ionization rate and njj is the hydrogen density. The cooling per 
unit volume is proportional to the collision rate, or n^. The kinetic 
temperature, usually increasing with the ratio of heating to cooling, 
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Figure 1. The log of the computed gas kinetic temperature is shown as a 
function of the hydrogen density and (top) the cosmic ray rate relative to 
the galactic background and (bottom) the extra heating rate. 



is then proportional to rcR/«H- Lines of constant temperature tend 
to run at a roughly 45 degree angle corresponding to a constant 
r CR./ n H ratio. The curves turn over at the highest densities when 
collisional cooling is suppressed. 

In the "extra-heating" case the heating per unit volume, 
Gextra, is assumed to be independent of density for simplicity. The 
temperature-determining ratio of heating to cooling is then pro- 
portional to Gextra/^H- Contours of constant temperature run at a 
steeper angle as a result. 

In both cases the temperatures in the lower right hand corner 
of Figure Q] are close to the CMB. Temperature increases as the 
heating increases, reaching the highest values in the upper left cor- 
ner where the gas is hot and highly ionized. 

Figure [2] shows the log of the hydrogen molecular fraction. 



log molecular hydrogen fraction 
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Figure 2. The log of the hydrogen molecular fraction for the cosmic ray 
(top) and "extra-heating" (bottom) cases. 



The lower right corner is predominantly molecular. Gas becomes 
increasingly dissociated and eventually ionized along a diagonal 
extending to the upper left. In the cosmic ray case H2 predomi- 
nantly dissociates following suprathermal excitation to molecular 
triplet states. Excitations to the singlets dissociate roughly 10% of 
the time with the remai nder decaying into various levels within th e 
ground electronic state (Sternber g et all 19871 ; iDalgarno et alll999h . 
Thus cosmic rays dissociate and excite H2 even when the kinetic 
temperature is low. This is in contrast with the "extra-heating" case, 
where little excitation or dissociation is produced until the kinetic 
temperature becomes large enough for thermal collisions to drive 
the processes. 

Figure|3]shows the log of the emissivity Alt j (erg cm~ 3 s _1 ) of 
the 1-0 S(l) 2.121 ^.m line, one of the strongest H2 lines in the NIR 
spectrum. The "extra-heating" case is shown to illustrate the effects 
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Figure 3. The emissivity of the H2 2.121 umline for the "extra-heating" 
case. Contours are drawn at 1 dex intervals, ranging from a minimum of 
logAnj of -29 up to -20. 

of changing temperature. The emissivity is strongly peaked towards 
a narrow band running along lines of roughly constant temperature 
and H2 fraction. Neglecting second-order effects, the emissivity of 
the line, which is predicted to be collisionally excited, will be pro- 
portional to the product of the density n(H-2) and a Boltzmann factor 

4^joc„(H 2 )exp(-^/r) (1) 

where % ~ 7000 K is the excitation energy of the upper level of the 
transition. The gas is mostly molecular in the lower right corner of 
the diagram but the temperature is too low to excite the transition. 
The temperature is high in the upper left corner but there is little 
H2. The line is emitted efficiently across a narrow band where the 
product of the H2 density and Boltzmann factor is large. The ridge 
of peak emissivity occurs at a gas kinetic temperature of roughly 
2000 K. 

This result can be generalized. Because of the form of equa- 
tion Q] a Boltzmann factor that increases exponentially as temper- 
ature rises, and an H2 density that decreases as the temperature 
increases, the emissivity of a collisionally-excited line will peak 
when the kinetic temperature is roughly half of the excitation en- 
ergy. The large changes in emissivity shown in Figure [3] suggest 
that, if gas exists across the parameter space, a linear detector like 
a spectrometer will pick out only the peaks. Simila r effects are 
prese nt in emission lines of active galactic nuclei dBaldwin et all 

GUI). 

We now concentrate on the cosmic-ray cas e since their densi- 
ties h ave been measured in regions of the flow jSanders & Fabian 
120071) . Figure [4] shows the emissivities, the emissi on per unit vol- 
ume ( erg cm~ 3 s _1 ), of all H2 lines detected in I Johnstone et all 
J2007h . These would be multiplied by the cloud volume, or equiv- 
alently the mass of H 2, to obtain line luminosities. We adopt the 
Johnstone et al (2007) density of ;ih = 10 4 cm~ 3 and vary the 
cosmic-ray density. This could occur if gas were trapped in mag- 




3 4 5 

log cosmic ray / background 



Figure 4. The pred icted emissivities of all H2 lines observed by 
Johnstone et al ( 2007) are shown as a function of the cosmic-ray density. 
The density is «h = 10 4 cm~ 3 . The low-excitation rotational lines are (red) 
28.21 |am (solid), 17.03 [jm (long dash), and 12.28 |am (short dash). The 
mid-excitation (2000 K < E v < 8000 K) lines are (green) 2.121 (solid), 
2.033 (am (long dash), 2.223 yan (short dash), and 2.423 \xm (dash dot). 
High excitation {E v > 8000 K) lines are (blue) 1.748 yxm (solid), 1.891 (jm 
(long dash), and 1.957 \xm (short dash). The emissivity of each line peaks 
at a temperature that is proportional to its upper level energy, as found by 
Ijohnstone et~a3 J2007I) . 

netic field lines, so that «h is constant, but exposed to a range of 
cosmic ray fluxes, perhaps due to attenuation or shielding. This cor- 
respo nds to a vertical slice through Figures 1-3. 

Johnst one et all d2007l Tablell) demonstrate a correlation be- 
tween a level's energy and its population temperature 7p 0p , the tem- 
perature corresponding to the population distribution. If the lev- 
els are thermalized then T pop is the kinetic temperature. Figure [4] 
shows the predicted emissivities of all lines reported in that paper. 
If the levels are thermalized then the population temperature will 
be weighted towards r peak , the kinetic temperature where the emis- 
sivity peaks. 

The distribution of em issivities shown in Figure |4]reproduces 
the Ijohnstone et 

3 d2007h correlation. The observed (Jp op ) and 
predicted (7p eak ) peak temperatures are, for low-excitation levels, 
r pop = 300 ±20 K, 7p eak = 610 ±450 K, for intermediate excita- 
tion, r pop = 1730 ±250 K, r peak = 2300 ±400 K, and for high exci- 
tation, 7/ pop = 2580± 150 K, r peak = 2900±400 K. This agreement, 
generally within the scatter, is partially fortuitous because the pre- 
dicted kinetic temperature depends on the gas metallicity, which we 
simply set to galactic ISM abundances and depletions. The overall 
trend, where population temperature correlates with excitation tem- 
perature, is inescapable for gas with a non-radiative energy source. 



3 CONCLUSIONS 

The comparison between the relative intensities of H2 and H I lines 
in the Orion star-forming reg ion and the cooling-flow filament stud- 
ied by Ijohnstone et all J20071) (Table [B shows that the H2 lines are 
far stronger in the filament. Starlight fluorescence produces only 
weak H2 emission since only a small portion of the Balmer contin- 
uum can be converted into H2 lines. This strongly suggests that the 
lines are formed by more efficient, perhaps non-radiative, processes 
in cooling flows. 

Both cosmic rays and dissipative MHD waves are candidates 
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as extra energy sources. Powerful selection effects introduced by 
the microphysics cause each emission line to have a peak emissiv- 
ity at a temperature related to its upper level energy. Both the great 
strength of the H2 lines and the observed correlation between ex- 
citation potential and population temperature result. Actually, any 
heating source that simultaneously heats and dissociates H2 would 
produce similar effects. 

The selection effects cause a collisionally-excited H2 spec- 
trum to masquerade as classical starlight-induced fluorescence. In 
the radiative case the lowest-energy levels are collisionally popu- 
lated and their populations indicate the kinetic temperature. Higher 
levels are fluoresce nce excited and have a physically meaningless 
higher temperature ( Ivan Dishoeck 2004). A similar range of popu- 
lation temperatures is produced in the non-radiative cases consid- 
ered here, but all of the temperatures are real. 

The observed cosmic ray densities are large enough to pro- 
duce the observed emissi on if nn is lower than usually assumed. 
Sanders & FabiarJ d2007l) find a cosmic ray energy density of 
roughly 15 eV cm~ 3 for c entral regions of the Perseus cluster. For 
comparison Webber ( 1998) finds the density of relativistic electrons 
in the local interstellar medium to be 0.2 eV cm~ 3 . The ratio of 
cluster to local ISM CR ionization rates, the vertical axis in most of 
the plots in this paper, is of order 10 3 . The predictions scale roughly 
as the ratio rcR/«H- The peak emissivity for high-excitation lines 
(Figure[4]l, occurs at rrjR as 10 5 5 for a density of tin = 10 4 cm~ 3 . A 
densit y of wh = 10 cm~ 3 w ould have similar spectra at rcR ~ 10 3 . 
Since[fohnstone et al ( 2007) was published, new calculations show 
that some H? collision rates were underestimated by roughly 2 dex 
dWrathmall et~alll2007l) . Simple scaling relations suggest that this 
would lower «h an d the gas pressure by a similar amount. This 
would solve the long-standing puzzle posed by the high gas pres- 
sures associated with the high density d Johnstone et alll2007r) . The 
filaments' gas pressure would then be in line with the pressure in 
the surrounding hot plasma. 

The emissivities produced by the no n-radiative proce s ses w e 
discuss are similar to those assumed by John stone et all d2007h . 
Hence the mass of H2 required to account for the observed emis- 
sion will be similar. T he total mass in warm H2, roughly 10 Mq 
d Johnstone et all|2007|) , is a tiny fraction of the mass of cold gas, 
~ 4 x 10 10 M Q , reported bv lSalome et~a3d2006l) . This suggests that 
the warm H2 may be produced in small regions which are affected 
by the localized deposition of energy in the forms of dissipative 
MHD waves, cosmic rays, or shocks. We have shown here that such 
sources of energy can account for the observed spectra. 
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